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Abstract. We derive the formulae for the time variation of the gravitational 
"constant" G and of the fine structure "constant" a in various models with extra 
dimensions and analyze their consistency with the available observational data for 
distant supernovae. We find that the reported variation of a translates into a small 
variation of G that makes distant supernovae to appear brighter, in contradiction with 
recent observations of high z supernovae. The significance of these results within the 
framework of some cosmological scenarios is also discussed. We find, however, that 
the magnitude of the effect is not large enough to safely discard the models with extra 
dimensions studied here. 
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1. Introduction 



One of the most challenging issues of modern physics is the quantization of the 
gravitational interaction. Within this theoretical framework it is worth noticing that 
although the controversial issue of the variation of the fundamental couplings dates back 
to the large-number hypothesis of Dirac p], it has recently become a subject of intensive 
experimental and theoretical studies — see |2] for an excellent review. Modern theories, 
like the string/M-theory or brane models, do not necessarily require a variation of the 
fundamental constants but they provide a natural and self-consistent framework for such 
variations. Most of such theories assume the existence of additional dimensions of the 
space-time and, therefore, contain a built-in mechanism allowing for the time variation 
of the observed couplings in four dimensions. This feature is rather easy to understand. 
Within the multidimensional approach the interactions are described by a fundamental 
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theory formulated in 4 + d dimensions, and the conventional four-dimensional theory 
appears as a result of the dimensional reduction. Couplings in four dimensions are 
determined by a set of a few constants of the multidimensional theory and the size R 
of the space of extra dimensions. The multidimensional constants are assumed to be 
genuinely fundamental and, consequently, do not vary with time. On the other hand, 
in the cosmological or astrophysical context it is natural to assume that R varies with 
time, very much in the same way as the scale factor of the three-dimensional space 
does. This leads to the time variation of parameters of the effective four- dimensional 
theory, like the gravitational constant G and fine structure constant a [3]-[7j. Moreover, 
since their time dependences are given by the same scale factor R these variations may 
turn out to be correlated, and it is natural to wonder what kind of relation is predicted 
in multidimensional scenarios. The effect of these correlations for various couplings in 
the framework of theories of unification has been discussed in ;8J. It has also been 
suggested that these correlations may unveil the mechanism of the time variation of 
the couplings and, in this way, give insight into the underlying fundamental theory. It 
is important to emphasize that a variation of the fundamental couplings may lead to 
deeper consequences or effects in the theory of particle interactions. In particular, it 
has been shown that a variation of fundamental couplings can be intrinsically related 
to the Lorentz and CPT violation 

There have been several attempts to measure the rate of variation of the 
fundamental constants but many of them have yielded just upper bounds on the 
absolute value of the rate of change. However, very recently, the rate of variation 
of the fine structure constant has been measured using high resolution spectroscopy of 
QSO absorption systems. To be specific a non-zero detection of the variation of the 
fine structure constant, Aa/a ~ — 10~^ at z ~ 1.5, has been reported (TUl. Although 
this result is still the subject of strong debate it is worth studying its theoretical and 
observational consequences. In addition, the observations of high redshift {z > 0.1) 
Type la supernovae [HI [12] and the analysis of the CMB strongly suggest a flat, Qr = 0, 
Universe ^3], with a mass density Qm — 0.3 and a non- vanishing cosmological constant 
fiyv — 0.7. This result, in turn, motivated a considerable amount of papers looking 
for viable alternatives. Amongst these alternatives the variability of the gravitational 
constant was one of the proposed alternatives [Hj to reconcile the observational data of 
distant supernovae with an open f^A = Universe. 

As it has been argued in JH], models in which the time variation of couplings 
— and in particular of the fine structure constant — is generated by the dynamics 
of a cosmological scalar field face serious difficulties. Namely, using rather general 
arguments, these authors show that in order to explain the observations of [TUj such 
models require an extremely precise fine-tuning that cannot be explained by any known 
mechanism. The problem is essentially the huge back-reaction produced by varying 
couplings on the vacuum energy. This difficulty could be intimately related to the 
long-standing cosmological constant problem. Hence, its satisfactory solution could also 
provide a mechanism to suppress the enormous variation of the vacuum energy due to 
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the time variation of a, but this is beyond the scope of this paper. Having adopted this 
point of view, we do not discuss this issue furthermore in the present work. Instead, 
the aim of the our work is to derive the relations between G/G and a/a in various 
models in (4 + d) space-time dimensions and to confront them with the observational 
data, namely to compare the detected variation of the fine structure constant [TO] with 
the existing bounds on the variation of G at cosmological distances. The models we are 
going to consider include the classical Kaluza-Klein models — see |161 [T7j for reviews 
on the subject — models with multidimensional gravitational and Yang-Mills fields 
|18| IT^ . and the Randall-Sundrum model with two branes and gauge and fermionic 
fields propagating in the bulk |2ni^[22]- The paper is organized as follows. In Section 
12] formulae for the time variation of G and a in various models are derived. In Section 
01 we use our theoretical results and the experimental value for a/a ^H] to obtain an 
estimate of the rate G/G, which is then confronted with the data obtained from distant 
Type la Supernovae [TTl [T2] . Finally, in Section HI conclusions and some discussion of 
the results are presented. 

2. Time variation of couplings in multidimensional theories 

In this section we consider theories with extra spatial dimensions. For the sake of 
simplicity we assume that the geometry of the space of additional dimensions is described 
by just one scale factor R. In most of the cases the generalization to geometries with 
a few scales is straightforward and does not bring any qualitatively new features to 
the effect we are going to study. Moreover, the theories studied here are assumed to 
be part of a cosmological scenario with varying R. The specific form of the function 
R{t) depends on the details of the scenario — see, for instance, jSH]. Here we simply 
assume that the evolution of the multidimensional Universe at the cosmological scale is 
described by some scenario which predicts certain background metric with varying scale 
R{t). Another possible class of scenarios would be the ones postulating a first order 
phase transition in the Universe. However, we will not consider this last possibility here 
and we limit ourselves to multidimensional theories describing the time variation of G 
and a which fulfill the following conditions: (1) the four- dimensional effective theory — 
that is, the dimensionally reduced theory — includes both Einstein gravity and Maxwell 
electrodynamics; and (2) the time variation of the scale factor R{t) of the space of extra 
dimensions leads to the time variation of both the gravitational constant G and the fine 
structure constant a. To calculate G/G and a/ a we focus only on the gravitational 
and electromagnetic sectors of the reduced theory. These time variations as well as that 
of R{t) are supposed to be slow enough in comparison with the phenomena described 
by these sectors — namely, the electromagnetic processes at the microscopic level. Our 
analysis will be rather general and will not rely on any particular function R{t). 
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2.1. Kaluza- Klein theories 



In two pioneering papers j24j T. Kaluza and O. Klein formulated the essential elements 
of the multidimensional approach to the description of fundamental interactions which 
was later called the Kaluza-Klein approach. They considered the pure Einstein gravity 
in five-dimensional space time M"^ x described by the multidimensional metric 
tensor qmn and showed that the sector of zero modes of the dimensionally reduced 
theory includes the four- dimensional gravity and Maxwell theory. Here is the four- 
dimensional Minkowski space-time, and is the circle. Later this construction was 
generalized |2H] to more general compact spaces of extra dimensions. In these cases the 
reduced theory contains Einstein gravity and Yang-Mills fields with the gauge group 
determined by the isometry group of the space of extra dimensions. In the standard 
setting the multidimensional Kaluza-Klein theory is the pure Einstein theory on M4 x Kd 
with the action given by 

S= [d'^'x^g—^ 7^(^+'^). (1) 

Here M4 is a (curved) four-dimensional space-time, Kii is a compact manifold of extra 
dimensions, g = detgMN, (M, N = 0, 1, 2, . . . , 3 + c?), TZ^^+^^ is the scalar curvature in 
M4 X Kd, and ^(4+^) is the multidimensional gravitational constant, which is assumed 
to be truly constant and does not depend on time. According to the procedure of 
dimensional reduction, to obtain the four- dimensional effective theory, firstly the fiu- 
components (/i, = 0, 1, 2, 3) of the metric tensor are identified as the four- dimensional 
metric tensor. Certain combinations of the rest of the components, ^^m, gnu and gmn 
with m,n = 4, . . . ,3 + d are identified as gauge field multiplets and scalar fields 
(pmn- Secondly, the mode expansion of all these fields is performed — see, for example, 
|17j . The coefficients of the expansion depend only on x^ and are interpreted as four- 
dimensional fields. In general there is an infinite number of them but here we are 
interested in the sector containing only zero modes of the mode expansion. Its action is 
given by 

where G{t) = (7(4) (t) is the four-dimensional gravitational constant. The parameters 
gi(t) are the gauge couplings, and the index i labels the simple subgroups of the gauge 
group. The dimensional reduction of the initial Kaluza-Klein action S, given by Eq. 
(jT]), yields in addition to So, given by Eq. terms including non-zero modes of 

the gravitational, gauge and scalar fields as well as terms proportional to {R/Ry. 
The scalar fields usually give highly non-linear interaction terms and are coupled non- 
minimally to the gravitational and gauge fields. For the sake of simplicity, the scalar 
fields are supposed to be frozen out and their contribution is neglected. Identifying the 
gravitational and gauge couplings from the action Sq for the zero modes we obtain the 
following expressions for G{t) and gi{ty in terms of G'(4+d) and the radius R{t) of the 
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space of extra dimensions: 

Git) ^ (3) 

where V^(t) ~ R{tY is the volume of the space of extra dimensions and ki are coefficients 
which depend on the isometry group of K^. Equations (jH)) and (jD) should be regarded 
as leading order approximations. The terms omitted in the reduced action will give 
sub-leading corrections, in particular through loop effects. As it has been said above, 
we assume that the dimensionally reduced theory includes the electrodynamics. Then 
the fine structure constant a{t) is given by a linear combination of gf{t), the specific 
relation depending on the model, in particular on the gauge group and the scheme of 
the spontaneous symmetry breaking. From Eq. Q it follows that 

where ki is some constant. Therefore, the time variation of G and a is determined by 
the function R(t). Since Vd/Vd = d{R/R), we get 

-=-(^+2)-. (7) 

As a consequence the time variation of the fine structure and gravitational constants 
are related by 

a _d + 2G 
a d G 

2.2. Einstein- Yang-Mills theories 



(8) 



Consider a theory in the (4 + (i)-dimensional space-time M4 x that includes gravity 
and the Yang-Mills field with the action 



S = I d^+'^xJ-g 



1 7^(4+rf) + ^l_TrFMNF'''' 



(9) 



where, as above, G(4+d) is the multidimensional gravitational constant, and g{i+d) is 
the multidimensional gauge coupling. Both are supposed to be constant in time. The 
dimensionally reduced theory includes the Einstein gravity and the four- dimensional 
gauge fields with an action similar to that of Eq. and, in addition, scalar fields with 
a quartic potential. The explicit form of the dimensionally reduced theory depends on 
the topology and geometry of the space of extra dimensions and the multidimensional 
gauge group. The case of Kd being a homogeneous space was studied in detail in 
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the literature (TH] — see also for reviews on the subject. The four dimensional 
gravitational constant is given by 

Git) ^ S^, (,0, 

We assume that the gauge part of the initial multidimensional model is such that 
its dimensional reduction gives the bosonic sector of the electroweak Glashow-Salam- 
Weimberg model in four dimensions. Examples of this kind can be found in ^H] ■ Then 
the fine structure constant in the reduced theory is related to the multidimensional 
gauge coupling 51(4+^) by 

am ^ (11) 

where K2 is some constant factor. From these expressions the following relation between 
the time variations of G and a can be easily obtained: 

a Lr 

We would like to mention here that the same realtion appears in a ten-dimensional 
model obtained as a low energy limit of a string theory [U]. 

2.3. Randall- Sundrum model 

A different multidimensional setting motivated by string/M-theories was proposed and 
studied in [20] • The model is formulated in the five-dimensional space-time x Ki 
with the fifth dimension compactified to the orbifold Ki = S^/Z2 of radius R — see 
[SE] and |2Zj for reviews on the subject. In the initial version of the Randall-Sundrum 
model with two branes located at the fixed points of the orbifold, known also as the RSI 
model, only gravity propagates in the five-dimensional bulk. The background metric 
solution is given by 

ds^ = e-^^^'^^'^^-^^^^dx^'dx" + i?2rf02 ^13) 

Pm I2HI 121]; where is the coordinate of the orbifold (0 < < tt), r/^,^ is the four- 
dimensional Minkowski metric tensor, and A; > is a parameter of the dimension of mass, 
its value being fixed by the brane tensions. Our three-dimensional space is identified with 
the brane with negative tension at = vr. The fields of the Standard Model are assumed 
to be localized at this brane. Due to the warp factor in the background metric, Eq. (|T^. 
the RSI model provides an elegant geometrical solution to the hierarchy problem |20j . 
In addition it opens new phenomenological possibilities like the observation of higher- 
dimensional gravity effects in the current or future collider experiments PU] . 

The reduction formula expressing the four-dimensional Planck mass in terms of 
the fundamental (five-dimensional) mass scale M = (167rG(5))^^/'^ ~ is derived j^H] 
using a metric given by Eq. which corresponds to the Galilean coordinates on the 
physical brane at = vr. One gets 

- (14) 
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(see also j2Z|)- To generate the correct hierarchy between the Planck scale and the 
TeV-scale the product kR must be kR ^ 11 12. As before, let us now consider 
this model as a part of a more general cosmological scenario with a slow varying scale 
factor R(t). The background solution (fT^ must be modified accordingly, in particular 
its four- dimensional part takes the form of the Robert son- Walker metric with the scale 
factor a{t) multiplied by the conformal warp factor including R{t). Models of this type 
have been extensively studied in the literature — see, for example, |3T|- It is natural 
to expect that the reduction formula in the first approximation remains the same as in 
Eq. ()14|) but with R = R{t). From this equation it readily follows that 

167rM3 e^'^-^W - 1 levrM^ ' ^ ^ 

The time variation of the Newton constant is then given by 

^ = -2vrfci?(t)|^— 1^ ^ -2nkRit)^ (16) 

However, since the fields of the Standard Model are localized on the brane and do not 
depend on R the RSI model does not contain a simple mechanism for the variation of 
the fine structure constant. To describe this effect one has to consider bulk gauge and, 
perhaps, fermionic fields. Such models have been studied in a number of papers, see 
for example j^H 1221 • assume that, similar to the case of the Einstein- Yang-Mills 
theories studied in §2.2, the bulk gauge fields yield the gauge (and perhaps the scalar) 
sector of the Standard Model on the brane. The electromagnetic ?7(l)-field appears in 
a usual way as a part of this sector after the spontaneous symmetry breaking. Similar 
to Eq. pi|) . the fine structure constant turns out to be related to the multidimensional 
gauge coupling g^^) by 

where is some constant. From Eqs. (|T6|l and ()17|) we obtain the following relation: 

^ = I -e-^'^-^W) i - (18) 

a 2nkR{t)G^ ' 2nkR{t) C ^ ' 

We would like to mention at this point that as it was observed in [7j the effect of the 
time variation of couplings in brane-world scenarios is closely related to the resolution 
of the hierarchy problem. 



2.4- General remarks 

Let us make a few remarks. Expressions (0), pi|) . and ()17|) are classical, or tree- level 
relations. They define the fine structure constant a{M^;t) at the scale M/j = R~^{t). 
We have found that a{Mji;t)/a{MR;t) = p{R/R), where p is a constant. In order 
to relate a{Mji;t) to the correspondig value a{fi;t) at some low energy scale /x, for 
example at the electroweak scale /i = Mz, one should take quantum corrections into 
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account by using the renormalization group formulas for running couplings. By standard 
considerations one obtains a relation of the form 
1 1 



+ A\n{fi/MR) 



where A is a constant of order one [S2I- Though, in principle, the second term also 
contributes to the time variation of a{n;t), in fact its variation is dominated by the 
first term, a~^(M/j;t). The reason is quite simple. Let us take the point of view that 
the scale fi does not vary with time. In other words, the variation of the dimensionless 
quantity fi/Mfi{t) is determined by R{t). Then by a straightforward time derivation of 
the previous expression one obtains: 

a[jj;t) K it L J 

The first term in this expression is the time variation calculated in §2.1, 2.2 and 2.3 
for three specific classes of models; the second term is of order 0{a). Hence, it is sub- 
dominant and will be neglected from now on. If the scale /i is defined in a different way 
and appears to be time dependent, the form of the second term in Eq. (fT^ may change, 
but it will remain to be sub-dominant. A similar analysis was presented in 1^ and [S]. 

In order to summarize, Eqs. (jH)), ()12|). and (fTHj) can be written in the following 
general form: 

^ = m)% (20) 

a G 

where 

for the Kaluza-Klein theories, 
(5[R) = ■{ 1 for the Einstein- Yang-Mills theories, (21) 



1 



for the Randall-Sundrum-type model. 



2TTkR{t) 

Note that a /a oc G/G, that the constant of proportionality is positive, and that /5 ~ 1 
in the case of the Kaluza-Klein and Einstein- Yang-Mills theories, and (3 ~ 10^^ in the 
case of the Randall-Sundrum model. We would like to emphasize that the result given 
in Eqs. and (PT|) is robust, since does not depend on details of the models and the 
specific form of R{t). 

If the dimensional reduction of the multidimensional models gives a four- 
dimensional model of unification, then an analysis similar to that carried out previously 
in subsections 2.1-2.3 gives the time variation of a single coupling constant ogut- In 
this case, Eq. ()19|1 for the coupling running is different and relates the three couplings 
of the Standard Model with acuT- The electromagnetic coupling at the low energy scale 
/i is calculated in a standard way. Making an analysis similar to that of jH] it can be 
shown that the time variation in the leading approximation is related to aGUT/«GUT 
by: 

a(/i;t) ^ ftGUT 
t) «GUT 
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where the proportionahty constant is of order unity. 

Finally, let us also mention that there is a class of models with branes and large extra 
dimensions, the Arkani-Hamed-Dimopoulos-Dvali (ADD) models, which also provide 
a solution to the hierarchy problem and predict new observable effects with massive 
gravitons jSSj- However, similar to the case of the RSI-model, the ADD-model in its 
standard version with only the gravitational field propagating in the bulk does not 
describe the time variation of the fine structure constant. 



3. Observational constraints on a/a and G/G 

The detailed analysis of three distant {z ~ 3.5) quasar absorption line data sets has 
provided for the first time direct evidence that the fine structure constant a was smaller 
in the past ^Hj- In particular, the detection of a variation in a: 

Aa a(z) — at) , ^ c: , ^ 

= ^ ^ = (-0.66 ±0.11) X 10-^ (22) 

a ao 

where ao = Ci{zo) is the present day value has been reported. It is of importance for our 
analysis to realize that Aa = a{z) — ao < 0. We also emphasize that, since the error 
bars are considerably smaller than the reported value, this is not just an upper bound 
but a direct measurement of the rate of change of the fine structure constant. Whether 
this detection is genuine or is affected by systematic errors is still a matter of debate 
[3 Ell- Using Eqs. (pUj) and (PT|) it turns out that, provided that the observational 
determination of Aa/a is correct, the variation of G/G in the models studied in the 
previous section is positive and we face a smaller value for G in the past. Given a 
typical age of the Universe tu ~ 14 Gyr and assuming a constant rate of change it is 
straightforward to derive from Eq. (j22j) the following estimate 

G/G ~ +10"^^ yr-i (23) 

for the Kaluza-Klein and the Einstein- Yang-Mills theories, whereas it is a factor of 10^ 
larger for the Randall- Sundrum scenario. 

Let us see if this estimate matches the currently available experimental bounds 
on the variation of G. Several constraints on the local rate of change of G using 
the observation of lunar occultations and eclipses, planetary and lunar radar-ranging 
measurements, the evolution of the Sun, gravitational lensing. Viking landers or data 
from the binary pulsar PSR 1913+16 have been obtained up to now. First, it is 
important to realize that all these measurements give just upper bounds on the rate 
of variation of G. Among these measurements, the last one provided for many years the 
most reliable upper bound [^3] 

-(1.10 ± 1.07) X 10-^Vr"^ < G/G < 0. 

However, the best upper bound has been obtained using helioseismological data jHEj: 

-1.6 X 10-^V"^ < G/G < 0. 
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Figure 1. The Hubble diagram of distant supernovae, assuming the preferred 
cosmological scenario of the SCP, for G/G = (solid line), G/G — +10~^^ yr""'^ 
(dashed line). For G/G = +10"^*^ yr~^ the results are indistinguishable from those 
of the case with constant G. The observational data are taken from JI]. The inset 
shows an enlarged view of the region around z 0.5. See the text for further details. 

Note that all these upper bounds are local. At cosmological distances the best upper 
bound for the rate of variation of G comes from the Hubble diagram of distant type 
la supernovae. By taking into account that distant Type la supernovae appear to be 
dimmer than local supernovae , the following upper bound on the variation of G 

was obtained |H7] : 

-IQ-" yr"^ < G/G < at z ~ 0.5. 

We would like to stress that all these upper bounds — regardless if they are local 
or obtained at moderately high redshifts — are negative and, consequently, positive 
values of G/G seem to be not allowed by the present astrophysical data. Consequently, 
these upper bounds are then at odds with the estimate given by Eq. (jSBI) which was 
obtained from the recent determination of a /a, Eq. ()22|) . within the multidimensional 
models studied in Sect. |2] In the following we will elaborate on this. It is as well 
worth mentioning at this point that the Oklo natural nuclear reactor [3H1 IBHI severely 
constrains a variation of a but, again, at low z. 
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Figure 2. Same as in Fig. 1, but assuming a flat, matter-dominated universe, 

inM,nA) = (1.0,0.0). 

Fig. n shows the Hubble diagram of distant Type la supernovae for three values 
of G/G and for the preferred scenario of the Supernova Cosmology Project, namely 
{Qr,Qm,^a) = (0.0,0.3,0.7). For the sake of simplicity we have assumed that G/G 
remains constant within this z interval. As it can be seen, the effect of a positive value 
of G/G is to make distant supernovae to appear brighter. This behavior is just the 
opposite to what it is observationally found. However, as it can be seen in Fig. ^ the 
effect of a varying G at a rate given by Eq. ()23|) . which is the value derived from the 
observed variation of a, is practically indistinguishable from that of G/G = 0. For the 
purpose of illustration we also show the curve with G/G = +10^^^ yr^^. Note also 
that the case G/G ~ +10"^^ yr~^ predicted in the Randall- Sundrum- type model — 
see Eq. fl21|l — is bracketed by these two curves. There are other scenarios which have 
been generally used to explain an accelerating universe without the need of invoking 
a non-vanishing cosmological constant by attributing this behavior to the presence of 
extra dimensions. To this regard in Fig.|21we show the Hubble diagram for a flat matter 
dominated universe {Qr, Qm, Qa) = (0.0, 1.0, 0.0). As it can be seen in this figure there 
is no possible way to reconcile the multidimensional models considered in Sect. EJwith 
the observations no matter which the value of G/G > is. 
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The conclusions above can be put into a quantitative form. Using the 2a confidence 
contours for z = 0.5 obtained from the fit to the Hubble diagram of Type la 
supernovae given in |22j we have calculated bounds to G/G. For the currently favored 
cosmological scenario {Qji,Qm,^a) = (0.0,0.3,0.7) and for the fiat matter dominated 
case {flR, flhi, ^a) = (0.0, 1.0, 0.0) we have obtained the following estimates respectively: 

nR = 0.0, nu = 0.3, = 0.7 -1.4 ■ 10-" yr-^ < G/G < +2.6 ■ 10~" yr-\ 
= 0.0, = 1.0, ^A = 0.0 -2.9 ■ 10-" yr-i < G/G < -0.3 • 10"" yr'^. 

Here, as before, we used the typical value tu = Hq^ = 14 Gyr and assumed a 
constant rate of change of G. As it has been already discussed above, the scenario 
{Qr, Qm, Qa) = (0.0, 1.0, 0.0) does not seem to allow a positive G/G. The cosmological 
scenario with {Qr,Qm,^a) = (0.0,0.3,0.7), which is the preferred scenario of the 
SCP, is among the allowed ones. It can be shown that for a fiat Universe positive 
values of G/G are allowed only if ^a >,0.15. Another interesting case is, for instance, 
{flR,QM,^A) = (0.5,0.5,0.0), for which we obtain 

-2.3 • 10"" yr-^ < G/G < +0.3 ■ 10"" yr^^ 

at the 2a confidence level. 

4. Conclusions and caveats 

We have derived the formulae for the time variation of the gravitational "constant" G 
and of the fine structure "constant" a for three classes of models with extra dimensions. 
We have found that such variations are related and we have derived the explicit relation 
— Eqs. fl20|) and (pT|) — which does not rely on the specific form of the time dependence 
of the scale factor of extra dimensions R{t). For the classes of models considered in ^ 
the relative sign between a/a and G/G turns out to be model independent. Then, using 
these expressions and the reported variation of a based on the available data obtained 
from distant QSOs [TOj, we have derived the estimate G/G +10"^^ yr"^- This value 
of the time variation of the gravitational constant makes distant supernovae to appear 
brighter, in contrast with observations. However, the effect is too small to safely discard 
the classes of models with extra dimensions considered here. In fact, a positive rate of 
variation G/G ~ +1 ■ 10"^^^^ yr"^ has been predicted within a iV = 1 ten-dimensional 
supergravity and with non-dynamical dilaton 0. However, and as it has been shown 
in jH], when the dilaton dynamics are taken into account the rate of change of a is too 
small to be observed. Let us mention that in this case the relation between a/ a and 
G/G appears to be given by Eq. (fT^ . 

We have also computed the Hubble diagrams of distant Type la supernovae in the 
framework of these models for several typical cosmological scenarios and analyzed their 
consistency with the available observational data. We have found that, if a flat Universe 
is assumed, models with extra dimensions and positive values of G/G, in accordance 
with Eq. (j22I), can reproduce the observational data only if f2A^0.15. One should 
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however keep in mind that provided that the value of a/a derived from QSOs turns 
out to be a genuine detection and given that the corresponding estimate of G/G is very 
small, at the redshifts of interest {z ~ 1.5) the deviation with respect to the Hubble 
diagram for G/G = is also very small so the observations still leave much room for the 
multidimensional models studied here provided that Q\ 0.15. Conversely, for the class 
of models studied here to be able to reproduce the observational data with f^A <,0.15 
negative values of G/G are needed (see Fig. 2). However, the absolute value of G/G 
needed to fit the observations should be much larger than that obtained here from 
distant QSOs. 

A natural question which arises is whether the robustness of the observational 
results on a/a and G/G allows to safely discard the multidimensional approach for 
scenarios involving a flat Universe with Q\<^0.15, attending exclusively to the relative 
signs of a/ a and G/G. There are a few possibilities. One of them is that the 
observational data is not precise enough at the moment and does not allow to draw 
any definite conclusions. This is perhaps the simplest and most obvious explanation. 
The situation may improve when, for example, better or additional observational 
determinations of a/a are obtained or when experimental data from distant supernovae 
become available in the future (with missions like SNAP) to give more accurate bounds 
on G/G. Although this is indeed the most straightforward explanation, let us consider 
other possibilities. Firstly, it may turn out that our consideration of the time variation 
of the fundamental couplings is too rough and some effects which modify Eq. (j^Ull 
are missing. And secondly, it could be that the multidimensional models considered 
here may also be too simple and phenomenologically unsatisfactory and should be 
substituted by more elaborated ones. In any case, questioning the applicability of the 
multidimensional approach for the description of the fundamental interactions should 
be preceded by additional more detailed theoretical studies and should be confronted 
with more accurate observational data. 
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